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Heisenberg interactions are ubiquitous in mag¬ 
netic materials and have been prevailing in mod¬ 
eling and designing quantum magnets. Bond- 
directional interactions^El offer a novel alternative 
to Heisenberg exchange and provide the build¬ 
ing blocks of the Kitaev modeP, which has a 
quantum spin liquid (QSL) as its exact ground 
state. Honeycomb iridates, A2lr03 (A=Na,Li), 
offer potential realizations of the Kitaev model, 
and their reported magnetic behaviors may be 
interpreted within the Kitaev framework. How¬ 
ever, the extent of their relevance to the Kitaev 
model remains unclear, as evidence for bond- 
directional interactions remains indirect or con¬ 
jectural. Here, we present direct evidence for 
dominant bond-directional interactions in antifer¬ 
romagnetic Na2lr03 and show that they lead to 
strong magnetic frustration. Diffuse magnetic x- 
ray scattering reveals broken spin-rotational sym¬ 
metry even above Tn, with the three spin com¬ 
ponents exhibiting nano-scale correlations along 
distinct crystallographic directions. This spin- 
space and real-space entanglement directly man¬ 
ifests the bond-directional interactions, provides 
the missing link to Kitaev physics in honeycomb 
iridates, and establishes a new design strategy to¬ 
ward frustrated magnetism. 

Iridium (IV) ions with pseudospin- 1/2 moments form in 
Na2lr03 a quasi-two-dimensional ( 2 D) honeycomb net¬ 
work, which is sandwiched between two layers of oxy¬ 
gen ions that frame edge-shared octahedra around the 
magnetic ions and mediate superexchange interactions 
between neighboring pseudospins (Fig. la). Owin g to 
the particular spin-orbital structure of the pseudospirP^, 
the isotropic part of the magnetic interaction is strongly 
suppressed in the 90 ° bonding geometry of the edge- 
shared octahedrsP^^, thereby allowing otherwise subdom¬ 
inant bond-dependent anisotropic interactions to play 
the main role and manifest themselves at the forefront 
of magnetism. This bonding geometry, common to many 


transition-metal oxides, in combination with the pseu¬ 
dospin that arises from strong spin-orbit coupling gives 
rise to an entirely new class of magnetism beyond the tra¬ 
ditional paradigm of Heisenberg magnets. On a honey¬ 
comb lattice, for instance, the leading anisotropic interac¬ 
tions take the form of the Kitaev modeP, which is a rare 
example of exactly solvable models with nontrivial prop¬ 
erties such as Majorana fermions and non-abelian statis¬ 
tics, and with potential links to quantum computing. 

Realization of the Kitaev model is now being 
intensive^ sought out in a growing number of 
materialP^, including 3 D extensions of the honeycomb 
Li2lr03, dubbed “hyper- honeycomb”l 3 and “harmonic- 
honeycomb’®, and 4 d transition-metal analogs such as 
RuClsP^ and Li2Rh02p^. Although most of these are 
known to magnetically order at low temperature, they 
exhibit a rich array of magnetic structures including zig- 
za^HHm^ spiraP ^, and other more complex non-coplanar 
structureP^LU that are predicted to occur in the vicinity 
of the Kitaev QSL phasd^^^^, which hosts many degen¬ 
erate ground states frustrated by three bond-directional 
Ising-type anisotropies. All of these magnetic orders are 
captured in an extended version of the Kitaev model writ¬ 
ten as 

H=Y^ [KS]S] + JSi ■ Sj + r(5f + sfs?)], (1) 

Sj)i 

which includes, in addition to the Kitaev term K, the 
Heisenberg exchange J that may be incompletely sup¬ 
pressed in the superexchange process and/or arise from 
a direct exchange proces^, and the symmetric off- 
diagonal exchange term F, which is symmetry-allowed 
even in the absence of lattice distortions^^HH] This 
“minimal” Hamiltonian couples pseudospins S (hereafter 
referred to as ‘spin’) only on nearest-neighbor bonds 
(ij), neglecting further-neighbor couplings, which may 
be non-negligible. The bond-directional nature of the 
K and F terms is reflected in the spin components 
[a 7^ /3 7^ 7 G {x,y,z)] which they couple for a given 
bond ( 7 G x—, bonds)(Fig. la). For example. 
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FIG. 1. Magnetic easy axis and temperature dependence of the zig-zag order, (a) Honeycomb layers of in 
the monoclinic Bravais lattice. Green, yellow, and blue planes show Ir206 plaquettes normal to local x, y, and z axes (black 
arrows), respectively, which point along three nearest neighbor Ir-0 bonds in an octahedron. Ir-Ir bonds are labeled following 
the plaquettes they belong to. Na^ is not shown for clarity. Blue arrows show the spins in the static zig-zag order propagating 
along the b direction. Spins are antiparallel between the layers (not shown), (b) Illustration of the scattering geometry. Shown 
in blue is the scattering plane defined by the incident (k^) and outgoing (k/) wavevector (red arrows). Green arrows show the 
x-ray polarizations. The azimuth T is dehned as the angle between the a axis and the scattering plane, (c) T-dependence of 
the magnetic Bragg peak (blue hlled circle) intensity at (0 1 3.5) measured in the a-n' channel. Black hexagon is the Brillouin 
zone of the honeycomb net. Red solid line shows the best htting result to the data with © = 44.3°. Green and blue lines shows 
the calculated T dependence for 0 = 40° and 50°, respectively, (d-f) H, K, and L scans of the magnetic Bragg peak at (0 1 
6.5) for selected temperatures, (g) Temperature dependence of the correlation lengths along the a, 5, and c axes from Gaussian 
htting to the scans. Error bars represent the standard deviation in the htting procedure. The solid lines are guides to eye. 


the K term couples only the spin component normal 
to the Ir206 plaquette containing the particular bond. 
Despite these extra terms that may account for hnite- 
temperature magnetic orders in the candidate materials, 
the fact that the Kitaev QSL phase ha s a h nite window 
of stability against these perturbationJ^SEH calls for in¬ 
vestigation of competing phases and a vigorous search for 
the Kitaev QSL phase. 

Although the notion of magnetic frustration induced 
by competing bond-directional interactions is compelling, 
it remains a theoretical construct without an existen¬ 
tial proof for such interactions in a real-world mate¬ 
rial. Moreover, theories for iridium compounds based on 
itinerant electrons suggest alternative picture^^^^^. In 
principle, measurement of the dynamical structure factor 
through inelastic neutron scattering (INS) or resonant in¬ 
elastic x-ray scattering (RIXS) provides the most direct 


access to the Hamiltonian describing the magnetic inter¬ 
actions. However, a fully momentum and energy resolved 
dynamical structure factor thus far remains elusive for 
any of the candidate materials; RIXS suffers from insuffi¬ 
cient energy resolutiorP^ and INS is currently limited by 
unavailability of large-volume single crystalJ^^. In this 
Letter, we take a new approach using diffuse magnetic 
x-ray scattering to provide direct evidence for predom¬ 
inant bond-directional interactions in Na2lr03 via the 
measurement of equal-time correlations of spin compo¬ 
nents above the ordering temperature (Tn = 12—15 K, 
See Fig. SI). 

We start from establishing the spin orientation in the 
static zig-zag or dei^i^®^ below Tn, as shown in Fig. la, 
using standard resonant magnetic x-ray diffraction. In 
this measurement, the x-ray polarization projects out a 
certain spin component; the intensity depends on the spin 
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FIG. 2. Diffuse magnetic x-ray scattering intensities above Tn* (a) Intensity plots in the HK plane measured at T = 17 
K for selected azimuth angles summing tt-ct' and tt-tt' channels, sensitive to spin components parallel to ki and perpendicular 
to the scattering plane, respectively. For example, ^ = 0° measures sum of correlations Sxx and Syy. The dashed hexagon 
indicates the first Brillouin zone of the honeycomb net. (b) Spin-component-resolved equal-time correlations extracted from 
(a), (c) Spin-component-integrated equal-time correlations extracted from (a). Peaks are located at Q = ±(0 1), ±(0.5 0.5), 
and ±(0.5 -0.5). (d) T-dependence of the diffuse peak intensities for Sample^l (open symbol) and #2 (closed symbol). Solid 
lines show the calculated T-dependence for x-, y-, and z-zig-zag states shown in (e) for the tt-ct' and tt-tt' polarization channels 
summed, (e) Zig-zag orders propagating along three equivalent directions. Blue zig-zag is the static structure, and green and 
yellow zig-zags are generated by 120° rotation of the blue zig-zag. 


orientation via the relation I oc |kj -Sp for the ct-tt' chan¬ 
nel measured, where k/ is the scattered x-ray wave vector 
(Fig. lb). Figure Ic shows the intensity variation as the 
sample is rotated about the ordering wave vector Q = (0 
1 3.5) by an azimuth ang le which causes S to process 
around Q. Earlier studiefl^^ have established that S is 
constrained to lie in the ac plane, so this measurement 
of /(^) determines the spin orientation by resolving the 
tilting angle 0 of S with respect to the a axis. The best 
fitting result with 0=44.3° indicates that the magnetic 
easy axis is approximately half way in between the cubic 
X and y axes (Fig. la). This static spin orientation is a 
compromise among all anisotropic interactions present in 
the system, and is strongly tied to the magnetic struc¬ 
ture because of their bond-directional nature. To see this 
point, consider, for example, the K term: in the zig-zag 
structure propagating along b direction, where the spins 
are antiferromagnetically aligned on the 2 :-bond and fer- 


romagnetically aligned on the x-bond and ^-bond, a fer¬ 
romagnetic (antiferromagnetic) K favors spins pointing 
perpendicular to (along) the z-axis for a pair of spins on 
the z-bond, and along (perpendicular to) the x-axis and 
^-axis for the pairs on the x-bond and ^-bond, respec¬ 
tively. 

The zig-zag order is one of the many magnetic states 
(including the aforementioned spiral and non-coplanar 
structures) that are classically degenerate in the pure 
Kitaev limifP^ and comprise the micro-states in the QSL 
phase. Away from the pure Kitaev limit, depending on 
their energy separations, signatures of other magnetic 
states and their associated magnetic anisotropies may 
become observable in the paramagnetic phase through 
diffuse magnetic scattering. In particular, zig-zag orders 
propagating in two other directions, ±120° rotated from 
the static one, are expected for a honeycomb net with Cs 
symmetry. [The actual 3D crystal structure has only an 
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approximate C 3 symmetry because of a monoclinic dis¬ 
tortion, which singles out one propagation direction for 
the long-range ordered state (along b directio^out of 
the three possible under the ideal Cs symmetry^.] 

With other magnetic correlations possibly emerging 
at high temperature in mind, we follow the tempera¬ 
ture evolution of the zig-zag order. Figures Id-f show 
K, and L scans, respectively, of the magnetic Bragg 
peak at Q = (0 1 6.5) for selected temperatures. Fig¬ 
ure Ig shows the correlation lengths along the a, 6 , and 
c axes as a function of temperature. As the tempera¬ 
ture increases above Tn, the zig-zag correlations dimin¬ 
ish rather isotropically despite dominant 2D couplings 
in the honeycomb net. This 3D characteristic of the 
magnetic correlations contrasts with that of the quasi- 
2D Heisenberg antiferromagnet Sr 2 lr 04 , which displays 
2D long-range correlations well above T#il, and implies 
that the critical temperature in Na 2 lr 03 is limited by the 
anisotropic interactions rather than the interlayer cou¬ 
pling; the Mermin-Wagner theorem requires either the 
symmetry to be lower than SU( 2 ) or the dimension to 
be higher than 2D for a finite temperature phase transi¬ 
tion. The zig-zag correlations survive on the length scale 
of several nanometers (approximately 3 unit cells wide) 
above Tn, but the peak intensities drop by two orders of 
magnitude. To isolate such small signals from the back¬ 
ground, we used an experimental setup that maximizes 
the signal-to-noise ratio as described in the Methods sec¬ 
tion. 

Figure 2a maps the diffuse scattering intensity over a 
region in momentum space encompassing a full Brillouin 
zone of the honeycomb net, at several different T angles 
to resolve the spin-components (see Fig. S 2 ). These maps 
integrate the dynamic structure factor over the range 0 
<uj <100 meV, covering the entire range of magnetic ex¬ 
citations (see Fig. 4), and serve as an excellent approx¬ 
imation for the equal-time correlation Saa = {SqS^q) 
{a = x^y^z). When averaged over the three spin com¬ 
ponents, the intensity map (Fig. 2c) indeed shows three 
zig-zag correlations above Tn with peaks at Q = ±(0 
1), ±(0.5 0.5), and ±(0.5 -0.5) of equal intensities, con¬ 
firming the near-ideal symmetry. However, the spin- 
component resolved maps, shown in Fig. 2b, manifestly 
break the symmetry. The system is left invariant only 
when (^3 rotation is performed simultaneously in the real 
space and in the spin space, i.e. cyclic permutation of 
spin indices. This ‘global’ symmetry implies a strong 
entanglement between the real space and the spin space. 
Specifically, the full azimuth dependence of each zig-zag 
state, shown in Fig. 2d, follows well the curves simulated 
for spin orientation fixed relative to the propagation di¬ 
rection, as depicted in Fig. 2e. In other words, specifying 
a spin component amounts to fixing the momentum di¬ 
rection and vice versa. This one-to-one correspondence 
between the spin space and the real space is a direct con¬ 
sequence of the bond-dependent nature of the anisotropic 
exchange terms. 

Qualitatively, it is immediately seen that the 
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FIG. 3. Simulation of diffuse scattering using exact di- 
agonalization. The Kitaev-Heisenberg model including up 
to third nearest-neighbor Heisenberg interactions was consid¬ 
ered. interpolates between the pure Heisenberg model and 
the pure Kitaev model via Ji=J 2 =J 3=1 — ^ and A = — A 
ferromagnetic K with finite Ji,J 2 , and J 3 stabilizes the zig¬ 
zag state for most values of The black curve shows the 
spectral asymmetry defined as the ratio of spectral weight at 
Q =(0.5,0.5) to that at Q = (0 1). Images show equal-time 
correlations {SqSIq) obtained by exact diagonalization using 
a 24-cite cluster and plotted in the extended Brillouin zone 
for selected The correlations for y and z components (not 
shown) can be generated by ± 120 ° rotations of the images 
shown. 


anisotropic interactions dominate over the isotropic in¬ 
teractions and the system is very far away from the pure 
Heisenberg limit, in which case the spatial correlations 
must be spin-component independent with three zig-zag 
peaks having equal intensities by symmetry (as in the 
spin-averaged correlation shown in Fig. 2c preserving C 3 
symmetry). A measure of how close the system is to 
either the Heisenberg or the Kitaev limits is provided 
by the intensity ratio of the weakest peak to the two 
bright ones in the spin-component-resolved correlations 
(Fig. 2b). To quantify this measure, represented by a 
variable linearly interpolating between these two limits, 
requires specifying the minimal Hamiltonian, which 
is not precisely known. For an estimation at a semi- 
quantitative level, we adopt a simple Hamiltonian that 
neglects all anisotropic terms beyond the K term. (This 
in turn requires including further neighbor Heisenberg 
couplings ±2 and J 3 to stabilize the zig-zag ordeJ^, which 
we take to be equal to Ji for simplicity.) Figure 3 shows 
the simulated patterns for selected It is clear that the 
observed diffuse pattern is consistent with the simulated 
one for the large ^ limit. In fact, the observed intensity 
ratio of ^0.2 is even smaller than calculated (Fig. 3) for 
the largest ^ in the zig-zag phase, which confirms the 
predominant anisotropic interactions. 

Interpreted within this model, our calculations would 
imply that the system is very close to the Kitaev limit. 
However, it is becoming increasingly evident that other 
aniso tropic terms beyond the Kitaev interaction do play 
a role^^^. This is, in fact, evident from the static spin 
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FIG. 4. Resonant inelastic x-ray scattering spectra be¬ 
low Tn* RIXS spectra recorded at T=9 K and 4^=180°. Q 
= (0 1), (0.5 0.5), and (0.5 -0.5), shown as blue, green, and 
yellow filled symbols, respectively, marked on the Brillouin 
zone of the honeycomb net and color-coded with the spectra. 
At this ^ angle, S lies approximately along k/ and tt-ct’ and 
tt-tt’ channels measure the two spin components transverse to 
S. 


not pointing along one of the cubic axes favored by the K 
term; all other anisotropic terms conspire to rotate the 
spin away from the principal axes. This in turn suggests 
that the zig-zag structure is further stabilized by other 
anisotropic terms. The zig-zag correlations survive at 
least up to ^70 K (see Fig. S3), which is in accord with 
the observation that coherent spin-waveJ^ disperse up 
to «5 meV. This energy scale coincides with the temper¬ 
ature scale («100 K) below which the magnetic suscep¬ 
tibility deviates from the Curie-Weiss behavioi®. This 
energy scale is, however, still far too small in compari¬ 
son to the energy («100 meV) spanned by the magnetic 
excitations (Fig. 4), suggesting that the zig-zag order is 


an emergent phenomenon. Despite that the macroscopic 
degeneracy in the Kitaev QSL phase is reduced down to 
three zig-zags, the high-energy Kitaev interactions leave 
their signature in the low-energy sector: the three spin 
components, each carrying its own zig-zag, compete and 
melt the long-range order at a temperature much lower 
than that suggested by the Weiss tem per ature (©w)^ 
leading to a large frustration parameter^ (= ©w/Tn) 
of ^8. 

Beneath a rather mundane static collinear antiferro¬ 
magnetic structure, the spins on the honeycomb lattice 
undergo fluctuations of a peculiar type that suppress the 
tendency to order. To reveal their full dynamics remains 
a challenge, but RIXS allows a glimpse into their uncon¬ 
ventional nature. Figure 4 shows the excitation spectra 
recorded in the ordered phase, in the scattering geome¬ 
try that probes two spin components transverse to the 
static moment. The soft excitations (cj^^O) are observed 
at Q = ±(0.5 0.5) and Q = ±(0.5 -0.5), away from the 
Bragg peak at Q = ±(0 1). This is a notable exception to 
the universality held in conventional magnets that spin 
waves emanate from Bragg peaks by virtue of the Gold- 
stone theorem, and magnetic anisotropy is manifested as 
a spin-wave gap, even in systems with extremely large 
magnetic anisotrop}!^. By contrast, the spin gap in our 
system is small (unresolved in our spectra and estimated 
to be smaller than 2 meV from INS datsP^ compared 
to the overall energy scale of the system, despite the 
fact that the magnetism is dominated by the anisotropic 
terms. Rather, the anisotropy is manifested as the sepa¬ 
ration of the long-wavelength spin waves from the Bragg 
peaks, which is a natural consequence of each spin com¬ 
ponent displaying its own real-space correlations. Our re¬ 
sults directly reveal the key building blocks of the Kitaev 
model in Na 2 lr 03 , and establishes a new design strategy 
for the long-sought quantum spin liquids via the bond- 
directional magnetic couplings. 
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Supplementary Information 
A. Methods 

A.l Single crystal growth characterization 

Single crystals of Na 2 lr 03 were grown following two 
different recipes using Na 2 C 03 flux (Sample 7 ^ 1 ) and 
self-flux (Sample # 2 ). For Sample 7 ^ 1 , a mixture of 
Na 2 C 03 and Ir 02 with a molar ratio of 50 : 1 was melted 
at 1050° C for 6 hours followed by fast cooling at a rate 
of 100 ° C/hour down to 1000 ° C, slow cooling at a rate of 
1 ° C/hour down to 800° C and furnace-cooling to room 
temperature in a sequence. Hexagonal pillar-shaped 
crystals with typical dimension of 0.2 mm x 0.2 mm x 
0.4 mm were obtained after dissolving Na 2 C 03 flux in 
acetone and water. For Sample 7 ^ 2 , powders of Na 2 C 03 
were mixed with 10 - 20 % excess Ir 02 and were calcined 
at 700° C for 24 hours. Single crystals were grown on 
top of a powder matrix in a subsequent heating at 1050° 
C. Plate-like crystals with typical dimension of 5 mm x 
5 mm X 0.1 mm were physically extracted. 

The powder x-ray diffraction patterns of both Sample 
7^1 and 7^2 were consistent with the crystal structure in 





the C 2 /m space group as previous reported^^. Sample 
7^1 had a slightly lower Tn=12 K compared to Sample 
7^2 with Tn= 15 K, as measured by SQUID magnetome- 
try (Fig. SI a) and by resonant x-ray diffraction through 
the magnetic Bragg peaks at Q = (0 1 n-|-|) (n: integer) 
(Fig. Sib). Sample 7^1 was found to be of multi domains 
but had a superior crystallinity with 0 . 1 ° mosaicity (as 
compared to 0.5° mosaicity of Sample 7 ^ 2 ) (Fig. Sic), 
and thus was used for the resonant diffraction experiment 
(polarization analysis and measurement of the magnetic 
correlation length.) Sample 7^2 was found to be of a 
single domain and was used for the RIXS measurement. 
Both Sample 7^1 and 7^2 were used for the diffuse scat¬ 
tering measurement and gave identical results. 


A.2 Resonant x-ray scattering 

Incident x-ray was tuned to Ir L 3 edge (11.2145 keV). 
The resonant x-ray diffraction was carried out at 6 ID-B 
beamline of the Advanced Photon Source. The polar¬ 
ization analysis was performed in the vertical scatter¬ 
ing geometry using a pyrolytic graphite analyzer prob¬ 
ing a-n' channel. The RIXS was performed at ID20 of 
the European Synchrotron Radiation Facility. The total 
instrumental energy resolution of 24 meV was achieved 
with a monochromator and a diced spherical analyzer 
made from Si (844) and a position-sensitive area detector 
placed on a Rowland circle with 2 m radius. The diffuse 
magnetic scattering was performed using the RIXS spec¬ 
trometers at 9 ID, 27 ID, and 30 ID (MERIX) beamlines 
of the Advanced Photon Source where a monochromator 
of 90 meV bandwidth was used for an order of magnitude 
higher incident photon flux than that from the Si (844) 
monochromator. In these experiments, horizontal scat¬ 
tering geometry was used with the 7r-incident x-ray polar¬ 
ization measuring the sum of tt-ct' and tt-tt' channels. The 
in-plane momentum resolution of the RIXS spectrometer 
was ±0.048 A“^. The use of the RIXS spectrometers re¬ 
jecting all inelastically scattered x-rays outside of the 100 
meV energy window centered at the elastic line led to a 
significant improvement in the signal-to-noise ratio. A 
typical counting time of 2 hours was required for a map 
shown in Eig. 2. 




theta (degree) theta (degree) 


FIG. SI. Characterization of Sample 7^1 and Sample 

#2. (a) Temperature dependence of magnetic susceptibility. 
Dotted lines indicate Tn- Black arrows indicate that the data 
were measured while warming after zero-field cooling, (b) 
Temperature dependence of the magnetic Bragg peaks, (c) 
Sample mosaicity. 


B. Extraction of the spin-component resolved 
equal-time correlators 

The x-ray scattering intensity measured without using 
a polarization analyzer contains contributions from both 
TT-cr' and tt-tt' channels, probing spin components along 
ki and perpendicular to the horizontal scattering plane, 
respectively. In other words, two spin components per¬ 
pendicular to k/ are measured in the 90° horizontal scat¬ 
tering geometry used. For example, when T=0°, the ^ 
local cubic axis points approximately along kj (Fig. S 2 ), 
and thus the scattering intensity measures the corre¬ 
lation Sxx^Syy. Likewise, T=120°(T=240°) measures 
zz(^^yy~^^zz^' Thcu, Sxx") ^yyi uud cau bc ex¬ 
tracted by solving a set of linear equations. 















FIG. S2. Scattering configurations, (a) 0°, (b) ^ = 

120°, and (c) ^ = 240°. The azimuth T is the angle between 
the a axis and the scattering plane. Thick black arrows de¬ 
note the local x, y, and 2 : axes in the IrOe octahedron, a axis 
is indicated by a thin black arrow. Cyan arrows are the in¬ 
cident and scattered x-rays with the wave vectors and k/, 
respectively. Green arrows indicate the x-ray polarizations. 
Q (dark blue arrow) is the momentum transfer. 


the diffuse map recorded at T=17 K (shown in Fig. 2a) 
apart from thermal broadenings. 



H(r. I. u.) 


b ^ = 0 ° 

50 K 




C. Temperature dependence of the diffuse mag¬ 
netic peak 

The short-range zig-zag order is observable at least up 
to T?^70 K (Fig. S3a). The magnetic correlation length 
(1.6-1.8 nm) along the a axis does not vary significantly 
in the measured temperature region. Fig. S3b plots the 
diffuse map at ^=0° at T=50 K, which is similar to 


FIG. S3. Diffuse magnetic x-ray scattering intensities 
above Tn* (a) H profiles of the diffuse peak at Q = (0 1). 
Shown in the inset is the magnetic correlation length along 
the a axis derived from Gaussian fitting (solid curves)of the 
data for Sample #1. (b) Diffuse map at T=0° at T=50 K. 
The dashed hexagon indicates the first Brillouin zone of the 
honeycomb net. 













